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Abstract 

A Bayesian analysis of the world's p(-f, K + )A data is presented. We adopt a Regge- plus- 
resonance framework featuring consistent interactions for nucleon resonances up to spin 
J = 5/2. The power of the momentum dependence of the consistent interaction structure 
rises with the spin of the resonance. This leads to unphysical structures in the energy 
dependence of the computed cross sections when the short-distance physics is cut off with 
standard hadronic form factors. A plausible, spin-dependent modification of the hadronic 
form factor is proposed which suppresses the unphysical artifacts. Next, we evaluate all 
possible combinations of 11 candidate resonances. The best model is selected from the 2048 
model variants by calculating the Bayesian evidence values against the world's p{"f,K + )k 
data. From the proposed selection of 11 resonances, we find that the following nucleon 
resonances have the highest probability of contributing to the reaction: Su (1535) , 5n (1650) , 
F 15 (1680), Pi 3 (1720), Di 3 (1900), Pi 3 (1900), Pn(1900), and F 15 (2000). 

1 Introduction 

The Regge-plus-resonance (RPR) framework provides a hybrid model that conjoins an isobar 
description of nucleon exchanges in the s channel and the exchange of Regge trajectories in 
the background channel. The members of a Regge trajectory share identical internal quantum 
numbers, such as strangeness and isospin, but have different total spins. More specifically, the 
spins and the squared masses of the Regge trajectory members are linearly related [1]. As 
the p(^f, K + )A reaction is dominated by background contributions at forward kaon scattering 
angles, only t channel Regge trajectories are introduced in the RPR model. The featuring 
Regge families of the RPR model are the if +(494) and iT* + (892) trajectories. By replacing the 
Feynman propagators of the K + (4Q4) and iv"* + (892) diagrams with the corresponding Regge 
propagators, the Reggeized background amplitude is obtained. 

In the s channel of the RPR model, individual nucleon resonances are exchanged. These 
nucleon resonances have half-integral spins and are modeled by effective Rarita-Schwinger (R-S) 
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fields.* The finite lifetime of a resonance is incorporated in the RPR model by replacing the 
real pole of the Feynman propagator, i.e. s — m R , with the imaginary pole s — m 2 R + imRTji. 
Here thr and Fr represent the mass and the decay width of the resonance, respectively. In 
order to ensure that the resonance contributions vanish in the high-s limit, a phenomenological 
hadronic form factor (HFF) is introduced. Such a form factor essentially "cuts off" the s channel 
amplitude beyond a certain energy scale. 

2 Consistent interactions: RPR-2007 versus RPR-2011 

In Ref. [2] a model was introduced that is now dubbed "RPR-2007" . In this model, only spin-1/2 
and spin-3/2 nucleon resonances are included. While spin-1/2 interactions are straightforward 
to deal with, spin-3/2 (in fact, spin J > 3/2 in general) couplings are more difficult to treat. 
The tensor-spinor representation of the R-S formalism inevitably invokes additional components 
for the R-S fields. These degrees of freedom are unphysical and are attributed to additional 
spin-1/2, 3/2, . . ., J — 1 components, next to the physical components of the spin- J R-S field. 
It is clear that these unphysical degrees of freedom should be eliminated from the transition 
amplitude in order to have reaction observables that are physically meaningful. 

In RPR-2007 so-called "contact-invariant" interaction Lagrangians are used to characterize 
the spin-3/2 couplings [2]. These Lagrangians contain one (two) "off-shell" parameter(s) for the 
strong (electromagnetic) vertex. There are hence three off-shell parameters for each amplitude 
that models a spin-3/2 resonance exchange. The HFF that is used in RPR-2007 is of the 
Gaussian form and has a common cutoff energy for all of the resonant amplitudes. The coupling 
constants, off-shell parameters, and the cutoff energy of the RPR-2007 model were optimized 
against the available, forward-angle experimental data at that time. The following resonance 
content was found for RPR-2007: Sn(1650), Pn(1710), Pi 3 (1720), Pi 3 (1900) and the missing 
resonance Z?i3(1900). Fig. 1 shows the RPR-2007 prediction for the p(^,K + )A differential 
cross section at three different photon energies in the lab frame. At forward kaon angles RPR- 
2007 is fully consistent with the data. At backward kaon angles, however, the model deviates 
significantly from the data. Moreover, an artificial bump is present and the situation worsens 
with increasing photon lab energy or, equivalently, increasing s. The observed elevation between 
the data and the RPR-2007 predictions can, to minor extent, be attributed to the fact that 
only cos 6^ > 0.35 data were considered when optimizing the RPR-2007 parameters. The main 
"The spin-1/2 R-S field is the Dirac field. Only fields with a spin J > 3/2 will be referred to as R-S fields. 
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Figure 1: RPR-2007 prediction and experimental data for the p(~/, K + )A differential cross section at photon 
lab energies of 1575, 1875, and 2175 MeV. The RPR-2007 model is optimized against the cost^ > 0.35 data, 
which is indicated by the arrows. More details about the calculations and the data can be found in Ref. [3]. 

cause, however, of the unphysical bumps at cos6* K < is the adopted inconsistent interaction 
Lagrangians for the spin-3/2 vertices. As mentioned earlier, these Lagrangians contain one or 
two off-shell parameters. Now, these parameters allow for physical couplings to the unphysical 
components of the spin-3/2 R-S field. Hence, they give rise to unphysical structures in the 
predicted observables, which cannot be eliminated for any combination of the values of the 
off-shell parameters. What is more, the off-shell parameters are actually fitted to the data and, 
consequently, so are the artificial structures. Clearly this is a peculiar situation and it has to 
be put right. 

In Ref. [4], a formalism was developed in which the interaction of R-S fields with funda- 
mental fields can be described in a consistent way. The consistent interaction Lagrangians are 
invariant under the so-called "unconstrained R-S gauge" and it was proven that the unphysical 
components of the R-S field decouple from the transition amplitude for this type of inter- 
actions. Moreover, the gauge-invariant Lagrangians do not contain any additional parameters, 
apart from the usual coupling constants. In Ref. [4] a novel, spin-dependent HFF was developed 
as well, dubbed the "multidipole-Gauss" HFF. This HFF is able to regularize the amplitudes 
involving consistently interacting R-S fields, something that is not feasible with the standard 
Gaussian HFF. The combination of consistent high-spin interactions and the accompanying 
multidipole-Gauss HFF, constitutes the foundations of a new RPR model, i.e. RPR-2011. 
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Figure 2: (a) The evidence values for the 2048 RPR model variants (blue circles) in function of the number of 
nucleon resonance parameters. The most probable model for a fixed number of parameters is indicated with a red 
square. Top right inset: evidence ratios relative to RPR-2011 for the models with the highest corrected evidence. 
The color coding refers to Jeffreys' scale: barely worth mentioning (orange), significant (yellow), strong to very 
strong (green) and decisive (white). More details can be found in Ref. [5] (b) The relative resonance probabilities 
for each of the 11 considered resonances. More details can be found in Ref. [3]. 

3 Bayesian inference of the RPR-2011 resonance content 

When it comes to specifying the set of nucleon resonances that has the most important contri- 
bution to the ^(7, K + )A reaction, one faces a lack of consensus between the various analyses 
for the reaction at hand. This disagreement becomes clear by inspecting Table I of Ref. [5]. 
The difficulty in determining the resonance contributions lies in the observed dominance of 
background (i.e. t channel) reactions in the p("y, K + )A process. 

The challenge for the newly developed RPR model is to determine the most probable model 
variant M (read: the most probable set of resonances) given the p^j, K + )A data {d^} of the 
last decade. The standard x 2 distribution for the RPR model variant space is only capable of 
specifying the set of resonances that offers the best description of the data {(4}. However, the 
X 2 distribution does not punish for the expansion of the set of resonances, i.e. the addition of 
extra model parameters. In fact, the absolute minimum of the x 2 hypersurface is most likely 
relocated, in general, by increasing the dimension of the model parameter space. Therefore, a 
model with a "large" number of contributing resonances is not a probable model. 

Bayesian inference offers a quantitative way of selecting the most probable model amongst 
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its possible variants. The Bayesian evidence Z for a specific set of resonances is defined as the 
probability of the data {c4}, given the model variant M, and can be expressed as 



with otM being the model's parameters. The integrand of the right-hand side of Eq. (1) is 
the product of the likelihood function C{olm) = P({dk}\o<-M, M) and the prior distribution 
it(olm) = P(&m\M). The prior distribution for the model's parameters olm is chosen to be 
a uniform distribution. Now, consider the two model variants M A and Mb- The probabilities 
for both models are given by P({dk}\M A ) and P({<4}|Mb). By using Bayes' theorem, the 
probability ratio P({dk}\M A ) / P({c4}|Mb) can be calculated as 



Since there is no prior preference for any of the two model variants, P(Ma) = P{Mb) and the 
probability ratio reduces to the evidence ratio. The most probable model is therefore the model 
with the highest evidence. Due to the non-Gaussian and correlated nature of the systematic 
errors of {dk}, the Bayesian evidence is underestimated. In Ref. [5] an approximate expression 
can be found for the corresponding corrected evidence, i.e. Z'. 

For the Bayesian analysis of the RPR-2011 model, 11 candidate resonances up to spin- 
5/2 are considered. So next to the consistent couplings, RPR-2011 differs from RPR-2007 
in the inclusion of spin-5/2 resonances. The candidate resonances are: 5*11(1535), 5*11(1650), 
£>i 5 (1675), Fi 5 (1680), £>i 3 (1700), Pn(1710), and Pi 3 (1720) (established resonances), Pi 3 (1900) 
and Fi5(2000) (less-established resonances), and Z?i 3 (1900) and Pn(1900) ("missing" reso- 
nances). Now, each of the 2 11 = 2048 model variants were optimized against {d^} (6148 data 
points to date) and the corresponding Z' value was calculated. The result is shown in Fig. 2(a). 
The model variant with the highest evidence is indicated as "RPR-2011". The resonances that 
are not included in RPR-2011 arc £>i 5 (1675), £>i 3 (1700), and Pn(1710). The best RPR model 
variant, i.e. the one that includes all of the 11 candidate resonances, is hence not the most 
probable model. The individual probability ratio for each of the candidate resonances is shown 
in Fig. 2(b). A positive ratio indicates that it is more probable that the resonance contributes 
to the reaction p(-j,K + )A than that it does not. Vice versa, a negative ratio implies that the 
possibility that the resonance contributes to the reaction is not supported by the data. There 
are only 3 resonances that have a negative probability ratio, which are exactly those that are 
not part of the resonance set of RPR-2011. 
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Figure 3: RPR-2011 prediction and experimental data for the p(~f, K + )A differential cross section at photon 
lab energies of 1275, 1575, 1875, and 2175 MeV. More details about the calculations and the data can be found 
in Ref. [5]. 

In Fig. 3 the RPR-2011 prediction for the ^(7, K + )A differential cross section is shown at 
four different photon lab energies. It is seen that the bumps at cos 6^ < (see Fig. 1) are no 
longer present and that RPR-2011 offers a good description of the data for the whole range of 
kaon scattering angles. 
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